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, We present a detailed study of the Higgs sector within an extension of the Minimal Su- 

d ■ 



persymmetric Standard Model that includes one Complex Higgs Triplet (MSSM+ICHT). 
The model spectrum includes three singly charged Higgs bosons as well as three CP-even 
(or scalar) and two CP-odd (or pseudoscalar) neutral Higgs bosons. We present an approxi- 



, mated calculation of the one-loop radiative corrections to the neutral CP-even Higgs masses 

^ . {t^h") s-'^d tlis couplings H^Z^Z'^ {i — 1, 2, 3), which determine the magnitude of the 

Higgs-strahlung processes e+e^ Z'^H^ . Limits from LEP2 arc then considered, in order 
to obtain bounds on the neutral Higgs sector. Further, we also include the experimental 
limits from LEP2 on e+e^ H^H^ and those on BR(t bH^) from Tcvatron, to derive 
bounds on the mass of the two lightest charged Higgs bosons (H^ and H^). Concerning the 
latter, we find some cases, where rnj^± ~ 90 GeV, that are not excluded by any experimental 
bound, even for large values of tan /3, so that they should be looked for at the Large Hadron 
Collider (LHC). 
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I. INTRODUCTION 

The Higgs spectrum of many well motivated extensions of the Standard Model (SM) often 
include charged Higgs bosons whose detection at future colliders would constitute a clear evidence 



of a Higgs sector beyond that of the Standard Model (SM) jl]. In particular, the 2-Higgs-Doublet- 
Model (2HDM, hereafter, of Type H), in both its Supersymmetric (SUSY) and non-SUSY versions 
3], has been extensively studied as a prototype of a Higgs sector that includes one charged Higgs 



However, a definitive 



boson pair (H^), whose detection is expected to take place at the LHC 
test of the mechanism of Electro- Weak Symmetry Breaking (EWSB) will require further studies 
aiming at pinning down the underlying complete Higgs spectrum. In particular, probing the 
properties of charged Higgs bosons could help to find out whether they are indeed associated with 
a weakly-interacting theory, as in the case of the most popular SUSY extension of the SM, the 
so-called Minimal Supersymmetric Standard Model (MSSM) [4], or with a strongly-interacting 
scenario, like the ones discussed recently j^. Ultimately, while many analyses in this direction can 
be carried out at the LHC, it will be a future International Linear Collider (ILC) saying the definite 
word about which mechanism leads to mass generation. Notice that these tests should also allow 
one to probe the symmetries of the Higgs potential and to determine whether the charged Higgs 
bosons belong to a weak doublet or to some larger multiplet. Among the latter, in particular, 
Higgs triplets have been considered Q], mainly to search for possible manifestations of an explicit 
breaking of the custodial SU{2)c symmetry, which keeps Veltman's so-called 'rho parameter' close 
to one, i.e. p ~ 1. Motivations to discuss Higgs triplets can also be drawn from models of neutrino 
masses [3] as well as scenarios with extra spacial dimensions Q]. Though most of the work has been 
within non-SUSY models jsl, there have also been studies of SUSY scenarios with complex Higgs 
triplets, such as in [2], where some phenomenological aspects of the Higgs sector were explored. 



Subsequent work in this model has been done in 



Decays of charged Higgs bosons have been studied in the literature [12], including the radia- 
tive modes W^jjW^Z^ (isl, mostly within the context of the 2HDM or its MSSM incarnation 
(including into SUSY particles [l^), but also for the effective Lagrangian extension of the 2HDM 



defined in [15[] and rnore recently within an extension of the MSSM with one Complex Higgs Triplet 
(MSSM+ICHT) All these activities are particularly relevant especially in view of the fact 
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that charged Hipgs boson decays can be exploited to detet.ine key pa_ of the ...detly™ 

Higgs sector [17]. Charged Higgs boson production at hadron colliders was studied long ago 



and, more recently, systematic calculations of production processes at the upcoming LHC have 



been presented [19[, including some higher order effects in QCD and SUSY QCD [20[]. Current 
bounds on the mass of the charged Higgs bosons have been obtained at Tevatron, by studying the 
top decay 4 — > 6 , which already eliminates some regions of the 2HDM and MSSM parameter 
spaces 



2ll ]. whereas LEP2 gives a model independent bound of mff± > 80 GeV [22l. |23|. 
In this paper we present a detailed study of the spectrum and discuss the phenomenology of 
the Higgs sector of the MSSM+ICHT model, i.e., the scenario that includes one complex Higgs 
triplet in addition to the usual MSSM Higgs content, namely two Higgs doublets. Our main focus 
will eventually be on the production and decay phenomenology of the charged Higgs states of the 
model. 

This article is organized as follows. In section II, we discuss the Higgs sector of this model, in 
particular, we present the charged Higgs boson spectrum and the inclusion of an estimated calcula- 
tion of the one-loop radiative corrections for the CP-even neutral Higgs sector. In this section, we 
also present a study of the couplings H^Z^Z^, which are modified by radiative corrections. Then, 
in section HI, we derive the expressions for the vertex H^ff (where / and /' are generic fermions 
with cumulative electromagnetic charge ±1) and we calculate the decay t — > b in the framework 
of the MSSM+ICHT model, also presenting numerical results for the most relevant charged Higgs 
Branching Ratios (BR's)^. (A comparison with latest bounds from Tevatron Run2 is also given 
therein.) A discussion of the main production mechanism at the LHC is presented in section IV. 
LHC event rates are given in section V. Finally, we summarize and conclude in section VI. 

II. THE CHARGED HIGGS SPECTRUM IN A SUSY MODEL WITH AN 
ADDITIONAL COMPLEX HIGGS TRIPLET 



The SUSY model with two doublets and a complex Higgs triplet (MSSM+ICHT) of [10| is one 
of the simplest extensions of the MSSM that allows one to study phenomenological consequences 
of an explicit breaking of the custodial SU{2)c symmetry 10|, llJ|- In the reminder of this section, 
we recap its main theoretical features. 



^ In the framework of the MSSM+ICHT the three charged Higgs states are denoted by with the convention: 



4 



A. The Higgs potential of the model 



The MSSM+ICHT model includes two Higgs doublets and a complex Higgs triplet given by 



he -^t 



(1) 



The Higgs triplet, of zero hypercharge, is described in terms of a 2 x 2 matrix representation: is 
the complex neutral field and , denote the charged fields. The most general gauge invariant 
and renormalizable Superpotential that can be written for the Higgs Superfields $1,2 and S is given 
by: 



W = X^i- S$2 + /Wi^i • ^2 + M2Tr(S^) 



(2) 



where we have used the notation $1 • $2 = £06^1^2- The resulting scalar potential involving only 
the Higgs fields is thus written as 

V = Vsb + Vf + Vd , 

where Vsb denotes the most general soft-Supersymmetry breaking potential, which is given by 

VsB = mll^if + ml\^2f + mlTr{J^^'E) 

+ [AX^i ■ S$2 + Si/xi$i • $2 + B2H2Tr(T,'^) + h.c], 



(3) 



Vp is the SUSY potential from F-terms 
Vf 
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Vd is the SUSY potential from D-terms 

^2 



Vd 
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In turn, the full scalar potential can be split into its neutral and charged parts, i.e., V = V^harged + 

^^eutral 



101, 



111 



Besides the Supersymmetry-breaking mass terms, mf {i = 1, 2, 3), the potential depends on 
the parameters A, ^i, /i2, that appear in Eq. (2), as well as the trilinear and bilinear terms, A 
and respectively. For simplicity, we will assume that there is no CP violation in the Higgs 
sector and thus all the parameters and Vacuum Expectations Values (VEVs) are assumed to be 
real. In the charged Higgs sector with the basis of states ((/)^, 0^*, ^2 ^ ^2*) a 4 x 4 

squared-matrix mass (M^)jj, i = j = 1, ...4. For the CP-odd Higgs sector with the basis of states 
{-^Imcpi, ■^Im(/)2, -^Im^'^) one obtains a 3 x 3 squared-matrix mass {Mp)ij, i = j = 1, 2, 3. 
For the CP-even Higgs sector with the basis of states {^Recp^, ■^Re4>2, "-"^^ ^ 

3x3 squared-matrix mass (Mg),;.,-, i = j = 1, 2, 3. The explicit expression of the resulting Higgs 



potential is given in Refs. 



111. 



B. Parameters of the model and definition of scenarios 

We can combine the VEVs of the doublet Higgs fields through the relation V£) = vf + V2 and 
define tan/3 = V2/V1. Furthermore, the parameters vd, vt, rri^^ and are related as follows: 

. _ 2lJ!D_ 



which implies that the /^-parameter is different from 1 at the tree level, namely, 

=i + m\ . (6) 

The bound on R is obtained from the p parameter measurement, which presently lies in the range 
0.9993-1.0006, from the global fit reported in Refs. 0, 23|. Thus, one has R < 0.012 and vt < S 



GeV. We have taken into account this bound in our numerical analyses. 

Thus, the Higgs sector of this model depends of the following parameters: (i) the gauge- 
related parameters {g, g' , v, R, tan/3); (ii) the Yukawa couplings (A, pi, P2) and (iii) the soft 
Supersymmetry-breaking parameters {A, Bi, B2). The gauge-related parameters can be replaced 
by the quantities {Gp^ a, m\\r, p, tan (3). For the numerical analysis to be realized in the remainder 
of this paper, we must make sure that the following theoretical conditions of the MSSM+ICHT are 
satisfied: (a) the global stability condition of the potential; (b) the necessary condition for having 
a global minimum and (c) the positivity of the mass eigenvalues of the full spectrum of charged, 
pseudoscalar and scalar Higgs bosons [10,] . 
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The parameter space analyzed here for the MSSM+ICHT is the same that was considered 
before in the speciahzed hterature 1^, [3], namely, we consider characteristic values below their 
perturbative limits. Small and large values of tan/3 are both considered. Typical cases for A, Bi, 
B2, fj-i and H2 are used to define the following scenarios. 

Scenario A. It is defined by considering Bi = ni = 0, B2 = —A and ^2 = 100 GeV while for A 
we shall consider the values A = 0.1, 0.5, 1.0. In this scenario it happens that the additional Higgs 
triplet plays a significant role in EWSB. 

Scenario B. This scenario is defined by choosing: B2 = ^2 = 0, Bi = —A, while for A we shall 
consider again the values A = 0.1, 0.5, 1.0. Most results will take /ii = 200 GeV, though other 
values (such as /ii = 400, 700 GeV) will also be considered. Here, the effects of the additional 
Higgs triplet are smaller, hence the behaviour of the model is similar to that of the MSSM. 

C. One-loop radiative corrections to the CP-even Higgs bosons masses in the 

MSSM+ICHT 

In some cases, within the Scenarios A and B, we will show that a very light CP-even Higgs boson 
appears at the tree level, with a mass around 10 GeV, which can even be as small as 0(0. 1) GeV. 
However, it is known that, in the MSSM, the inclusion of radiative corrections from top and stop 
loops can alter significantly the (lightest) neutral CP-even Higgs mass. Thus, we can expect that 
similar effects will appear here and, furthermore, one also needs to consider in the MSSM+ICHT 
a possible large correction from Higgs-chargino loops, which could lift the corresponding Higgs 
mass above current experimental bounds. This means that one needs to include all such radiative 
corrections in order to avoid misleading conclusions. We are also interested in discussing the neutral 
Higgs bosons masses here because of their possible appearance in charged Higgs boson decays. As 
it will be shown later, this effect is important for large regions of the MSSM+ICHT parameter 
space. 

The radiative corrections to Supersymmetric Higgs boson masses can be evaluated using the 



effective potential technique 



2J], which at one-loop reads: 



Vi{Q) = Vo{Q) + m{Q), (7) 

1 Af ^ S 

^'^^^Q) = M^2SirM\log-^--), (8) 

where Vo{Q) is the tree-level potential evaluated with couplings renormalized at some scale Q, 
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Str denotes the conventional Supertrace and is the mass matrix for the CP-even sector. As 
discussed in I2J], the radiatively-corrected Higgs mass-squared matrix is given by the matrix of 
the second derivatives of Vi with respect to the Higgs fields, which is written as a function of their 
self-energies. In the MSSM, we know that the most important contributions to the Higgs self- 
energies at the one-loop level come from the diagrams with the top quark (and its scalar partner) 
circulating in the loop, due to the large top Yukawa coupling. However, for very large values of 
tan P, the bottom-sbottom contributions can become non-negligible. Therefore, for our settings, 
the dominant contributions to the Supertrace in the MSSM-I-ICHT are due to the top-stop and 
bottom-sbottom loops. Within this approximation, it happens that the squared-mass matrix of 
the CP-even Higgs bosons only gets corrected along its (1, 1) and (2,2) elements, given as follows: 

4 

o TTl- 

(AM|)i,i = —Xlmllog^, 



Stt^ " " ~ m 



o 777,- 

(AM|)2,2 = ^Xlmllog-^, (9) 

where Xt^b are the Yukawa couplings and the D-terms are omitted. In short, in the MSSM-I-ICHT, 
the radiative corrections to the Higgs boson masses must include the dominant contribution from 
the top-stop and bottom-sbottom systems. For this, it is enough to suitably modify the elements 
(AM|)i^i, (AM|)2,2 to the squared-mass matrix M| of the CP -even Higgs boson. Furthermore, 
as intimated already, we must also evaluate the contribution from the fermionic partner of the 
Higgs Superfields, which includes the Higgs-Higgsino triplets, because there is a potentially large 
effect emerging in the calculation of the squared-mass matrix of the CP-even Higgs bosons when 
the parameter A is large. Similarly to the top-stop and bottom-sbottom corrections, we estimate 

^1)3.3 



that the correction from the Higgs-Higgsino only modifies the element {Mo)s.3 



f2^ _ 3 ,2_2 ,_™J± 



(AM^)3,3 = (10) 

where A is the Yukawa coupling that appears in the Superpotential of the Higgs Superfields and 
- within our approximation - we take 'm^± and mjj± as the mass scales of the lightest charginos 
and charged Higgs bosons, respectively, i.e., ?7i^± ~ ^x"^ "^-ff* — m^±. D-terms are omitted, 
as well as possible effects from stop, sbottom and Higgsino mixing. Previous studies of Higgs mass 
bounds of this model were considered by J. R. Espinosa and M. Quiros j^], who concluded that 
the lightest Higgs boson of the model satisfy the bound. 



^ mz'^cos2(2/?) + l/2{\^v'^ /ml)sm^{2[3). (11) 
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Thus, for values of A that are consistent with perturbativity, wich then imphes a bound of the order 
mjjo ^ 155 GeV. Throughout this paper we take values of A that do not saturate this bound. A 
more complete calculation of the radiative corrections at one-loop level for this model is in progress 



26|. 



The main consequence of these radiative corrections is that the lightest CP-even Higgs mass can 
be enhanced at such levels that it makes it possible to pass current experimental bounds from LEP2. 
Besides, the radiative corrections affect mainly the neutral Higgs bosons sector, in particular the 
production of the neutral scalar Higgs in e~^e~ collisions, which is the Higgs-strahlung processes 
e~^e~ — > HfZ^, whose cross sections can be expressed in terms of the SM Higgs boson (herein 



denoted by (p'gj^j) production formula and the Higgs-Z'^Z'^ coupling, as follows 



R 



H°ZOZO 



^Mfzozo 



(12) 



where 5^o^o^o is the coupling HfZ^Z^ in the MSSM-I-ICHT and g'^o ^q^o is the SM coupling 



Z^Z^, which obey the relation 



i=l 



ZOzo — d^o^^^zozo- 



(13) 



In particular, for our model the factor -R^o^o^o is given by: 



(14) 



where V^^ denote the zj'-elements of the rotation matrix for the CP-even neutral sector, which 
relates the physical states and the real part of the fields (j)^, cl)2, in the following way: 



V{i Vfi, V, 



12 "13 



yS yS yS 
"21 ^^22 ^23 

[ vil vi vi ] 



HO 



(15) 



where the V/j are modified by the one-loop radiative corrections to the CP-even sector of our 
model. For our numerical analysis of the Higgs mass spectrum in the MSSM-I-ICHT we consider 



the experimenta 
Higgs state Hf" 



limits on the charged Higgs mass from LEP2 and apply it to the lightest charged 



22 



23]. The bounds on the neutral Higgs bosons H^, H2 are expressed in terms of 



the LEP2 bounds for i?^ 



H?ZOzo 



27l | . We will show that this excludes large regions of the parameter 



space of the MSSM-I-ICHT model. This is summarized in Tables IB lIVI Herein, we define as the 
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"marginal regions" those cases that ahnost pass LEP2 bounds on the neutral Higgs, i.e., when 
n^H°2 ~ GeV and/or i?^,, ^o^o not consistent with experimental bounds but for which 
we expect that the complete calculation of the one-loop radiative corrections to the mass of the 
neutral Higgs boson in question could enhance its mass, thereby allowing it to eventually pass said 
experimental limits. 



D. Higgs masses: numerical results 

Let us consider first Scenario A. Figures [H [2] and Eljll [5] and [6]) show the results for 
charged (neutral) Higgs bosons masses as a function of tan (3, in the range 1 < tan (3 < 100, 
for the cases A = 0.1, 0.5, 1.0, while taking A = 200, 300, 400 GeV, respectively. Throughout 
this paper we shall assume that the numerical values for stop and sbottom masses, taken at the 
electroweak scale, are degenerated. The above results for charged Higgs massed is based on the 
tree-level analysis. Similarly, the coming results for the pseudoscalar masses is also based on the 
tree-level formulae. However, the masses of the neutral CP-even Higgs bosons is based in the 
previous discussion of one-loop radiative corrections to the Higgs masses. For the stop, sbottom 
and chargino masses we take as input the value m = 1 TeV. In Figure [1] we present the charged 
Higgs boson masses for A = 0.1. We can see that the lightest charged Higgs boson has a mass 
Tnjj± ^ m\Y± , which is not below the theoretical limit that one obtains in the MSSM. Similarly, 
Figure [2] shows the charged Higgs boson masses for the case A = 0.5, and again we have that 
Tnjj± ^ m\Y± is possible but only for large tan f3. Furthermore, here it is possible for both Hf^ and 
to be lighter than the top quark. Figure [3] shows the charged Higgs boson masses for the case 
A = 1: now the lightest charged Higgs boson has a mass in the range 100 GeV < nT,^± < 200 GeV. 

Figure m shows the neutral Higgs spectrum for the case A = 0.1, and we notice the presence of a 
light CP-even Higgs boson with 11 GeV < irifjo < 50 GeV, especially for low values of tan (3 (< 5), 
that at first sight it would seem excluded by the LEP2 experimental limits. In fact, when one 
compares the results for -R^o^o^o obtained for this model , which measures the strength of the 
Higgs-strahlung process, with the LEP2 bounds [^l, which require it to be less than 0.01, we 
conclude that this scenario is indeed excluded, as it is summarized in our Table HI We assume that 
the lightest neutral Higgs boson decays predominantly into bb mode. Similarly, Figure [5] considers 
the case A = 0.5, and again we find 11 GeV < rrifjo < 50 GeV for 1 < tan /? < 100. However, we 
find that, for tan (3 < 77, -^^o^o^o is within the range allowed by LEP2. There is also a region 
where 111 GeV < rrifjo < 114 GeV, which we identify as marginal. Finally, Figure [6] corresponds 
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to the case A = 1.0, and we find that 14 GeV < m^jo < 89 GeV, for 15 < tan/3 and, although 
iZ^o^o^o < 0.01, again we find that this is a marginal region because 111 GeV < rrij^o < 114 GeV. 

As a lesson from these figures, for the case A = 0.5, we find that the LEP2 limit on the charged 
Higgs mass allows cases where 'rn^± ^ m^±, while the neutral Higgs bosons (chiefly Hf) satisfy 
the experimental limits of LEP2. However, the case A = 0.1 is not a favorable scenario, because 
Rh^z^'zo does not satisfy the experimental bounds. In contrast, for A = 1.0, the charged Higgs 
boson masses are significantly heavier. A complete list of bounds for all cases considered within 
Scenario A is shown in Table HI 

TABLE I: Analysis of ^^o^o^o consistent with LEP. Wc consider experimental limits allowed by LEP2 for 
charged and neutral Higgs bosons, for Scenario A with A ~ 200, 300, 400 GeV and /i2 = 100 GeV. 



A = 0.1 


tan /3 < 5 


m„± 81 GcV 
11 GeV < m^o < 50 GeV 
111 GeV < m„o < 118 GeV 


0.15 < iJ^o^o^o < 0.8 
^H«zozo < 


Excluded by ^^o^jo^o 


A = 0.5 


tan /3 < 77 


79.8 GeV < m„± < 118 GeV 

12 GeV < m^o < 50 GeV 
111 GeV < m„o < 114 GeV 


0.002 < Rlo^ozo < 0-2 

< ^H«Z"ZO 


Allowed by R^ozozo' 
but marginal for ^J^o^o^o 


A = 1 


15 < tan (3 


89 GeV <m„± < 187 GeV 
^1 

14 GeV < m^o < 89 GeV 
111 GeV < mrra < 114 GeV 


^lozozo < O-Ol 
0.9 < R^o^azo 


Allowed by R%ozOzO' 
but marginal for ^J^o^o^o 



For Scenario B, Figures [3 [8] and [Q tjlOl [TT] and [T2]) show the charged (neutral) Higgs bosons 
masses, as a function of tan/3 in the range 1 < tan/? < 100, and for the cases A = 0.1, 0.5, 1.0, 
taking A = 200, 300, 0.1 GeV, respectively. The lowest value {A = 0.1) is designed in order to get 
charged Higgs masses below the top mass. Let us comment first the results found for the charged 
Higgs mass in the case A = 0.1, that appear in Figure [71 We can see that the lightest charged 
Higgs boson has a mass above 300 GeV for A = 200, 300 GeV, while even for A = 0.1 GeV, it has 
a mass above m^r±, but it is still lighter than the top quark. Similarly, Figure [5] shows the charged 
Higgs boson masses for the case A = 0.5. We find that, for A = 200, 300 GeV, rn^± > 300 GeV, 
while, for A = 0.1 GeV, the mass is still in the range 100 GeV < rnjj± < mt- In turn. Figure [9] 
shows the charged Higgs boson masses for A = 1. Now, we have that the lightest charged Higgs 
boson is heavier than the top quark, even for ^ = 0.1 GeV. 

Let us now discuss the neutral Higgs spectrum. Figure [TO] shows the case A = 0.1 for A = 200, 
300 GeV, where one finds that 60 GeV < m^^o < 110 GeV, for 1 < tan/3 < 100, but the region 
allowed by iJ^o^Q^g corresponds to 10 < tan (3, while the parameter area corresponding to nifjo ~ 
110 GeV is of marginal type. Then, the case A = 0.1 GeV gives neutral Higgs masses within the 
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TABLE II: Analysis of -R^o^o^o consistent with LEP. We consider experimental limits by LEP2 for charged 
and neutral Higgs bosons, for Scenario B with A — 200, 300 GeV and fii = 200 GeV. 



A = 0.1 


10 < tan /3 < 100 


294 GeV < m„± < 532 GeV 

m„o ^ 110 GeV 
300 GeV < m„o < 538 GeV 


0-99 < '^lozozo 
^lozozo < 001 


Marginally allowed by R^o^a^o 


A = 0.5 


1 < tan /3 < 100 


300 GeV < m„± < 1200 GeV 
^1 

100 GeV < m^o < 107 GeV 
290 GeV < m„o < 1200 GeV 

"9. 


0.99 < Rl,^,^, 
^lozozo < 0-001 


Marginally allowed by R^o^o^o 


A = 1 


1 < tan /3 < 100 


340 GeV < m„± < 1690 GeV 
^1 

104 GeV < mjjo < 176 GeV 

252 GeV < m„o < 1700 GeV 

H9. 


0.99 < Rl„^,^, 
^lozozo < 0-001 


Allowed by R^o^ozo 



range 14 GeV < rrijjo < 50 GeV, rrifjo 110 GeV, for which one finds that -R^o^o^o < 0.01. 
In this case we have a marginally allowed region. Figure [TT] corresponds to the case A = 0.5, for 
A = 200, 300 GeV, and now we have 100 GeV < rriHO < 107 GeV for 1 < tan /? < 100. However, 
we find that 0.9 < -R^o^o^oi so that this region is marginal for nifjo. For the case A = 0.1 GeV 
we have 16 GeV < rrifjo < 50 GeV, nifjo ~ 107 GeV, but also -R^o^o^o < 0.01, as we can see in 
Table Hill We consider that this is a marginally allowed region by LEP2. Furthermore, Figure [12] 
includes the case A = 1, for ^ = 200, 300 GeV, and we find 104 GeV < m^^^o < 176 GeV, with 
0.9 < -R^o^o^o) again a marginal region allowed by LEP2, namely when rrij^o < 115 GeV and 
10 < tan/3. Finally, in the case A = OA GeV, we have 14 GeV < nij^o < 46 GeV, rrifjo ~ 104 
GeV, and R^o^ogn^ 0.01, as we can see in Table HIH We also identify this as a possible marginal 
region for LEP2 [27]. 



As a lesson from this second set of figures, we can state that, for the cases A = 0.1, 0.5, 1.0, 
which leave the range tan/3 > 25 and A = 0.1 as acceptable, we have 84 GeV < mjj± < 200 
GeV, while the neutral Higgs bosons lay in the mass range disallowed by the experimental limits of 
LEP2. However, in the cases with A = 200, 300 GeV, we obtain masses of the charged Higgs boson 
heavier than the top quark. A complete list of bounds for all cases considered within Scenario B 
is shown in Tables [III and [Till 

Finally, in order to consider possible variations with /zi and the behavior in the limit A ^ 0, we 
present the charged Higgs boson masses as function of tan /3, in Figure [T3l for A = 0.5, A = and for 
/ii = 200, 400, 700 GeV. We can see that for all these cases ~ 130 GeV. The corresponding 

results for the neutral spectrum are shown in Figure [TH and in Table IIVI which indicate the 
presence of two light neutral Higgs states, which lay in the range: 10 GeV < rrifjo < 50 GeV, 97 
GeV < rufjo < 107 GeV, for 2 < tan^ < 6. In this case we also find 0.007 < Rhozozo < 0.06, 
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TABLE III: Same analysis of Table [U but taking A = 0.1 GeV. 



A = 0.1 


12 < tan /3 < 100 


84 GeV < m„± < 95 GeV 
14 GeV < rrijjo < 50 GeV 
m„o ~ 110 GeV 


0-99 < ^lozozo 


Allowed for -R^o^ozo 
Marginal region for R%ozOz<> 


A = 0.5 


4 < tan/3 < 100 


121 GeV < m„± < 129 GeV 
16 GeV < nijjo < 50 GeV 
m„o ~ 107 GeV 


^'hOzozo < O-Ol 


Allowed for R^o^ozo 
Marginal region for R^a^ozo 


A = 1 


27 < tan /3 < 100 


197 GeV < m„± < 200 GeV 

14 GeV < rrijjo < 46 GeV 
103 GeV < m^o < 105 GeV 


iJ^o^jo^o < 0.01 
0-99<«?fO^o^o 


Allowed for -R^o^ozo 
Marginal region for R%ozOzo 



TABLE IV: Same analysis of Tabic HH but taking A = GeV, ^J.l = 200, 400, 700 GeV and A = 0.5. 



A = 0.5 


1 < tan /3 < 6 


121 GeV < m„± < 130 GeV 
10 GeV < nijjo < 50 GeV 
97 GcV< m„o < 113 GeV 

"9. 


0.006 < R^o^ozo < 0.2 
0-76 < Rlozozo 


Allowed for -R^o^jo^o 
Marginal region for R\ozOzo 



which then gives a small region of tan (3 that could be allowed by the LEP2 data. In fact, for 
tan/? > 6, practically all these scenarios are ruled out by LEP2, because rrijjo < 10 GeV 

III. THE VERTEX H^ff AND THE DECAY t ^ H+ b 

In the two previous papers [ill . [3], some of us have studied the bosonic interactions of the 
charged and neutral Higgs bosons, now we shall discuss the corresponding fermionic interactions 
and their implications for charged Higgs boson production through top quark decays and via 
Higgs-strahlung off-top quarks. 



A. The Higgs boson coupling to fermions in the MSSM+ICHT 

As in the MSSM, also in this model only the scalar doublets are coupled to the fermions, so 
that the Lagrangian of the Yukawa sector has the following expression: 

Cvuk = -Xu[uPLUct>2 - uPLdct)^] - \d[dPLd^\ - dPLUCP^] + h.c, (16) 

where the parameters X^^d are related to the fermion masses via 

K = , Ad = . (17j 
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FIG. 2: Same as in Figure [U but taking A = 0.5. 
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FIG. 3: Same as in Figure [TJ but taking A = 1. 




FIG. 4: Mass spectrum of the neutral Higgs bosons, taking /i2 = 100 GeV and A = 0.1, for: A = 200 GeV 
(left), A = 300 GeV (center), A = 400 GeV (right). The solid lines correspond to scalar and the dashed 
lines to pseudoscalar eigenstates. 
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FIG. 5: Same as in Figure IH but taking A — 0.5. 
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FIG. 7: Mass spectrum of the charged Higgs bosons, taking iii — 200 GeV and A = 0.1, for: A = 200 GcV 
(left), A = 300 GeV (center), A ^ OA GeV (right). 
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FIG. 8: Same as in Figure [71 but taking A = 0.5. 
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FIG. 9: Same as in Figure [71 but taking A = 1. 
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FIG. 10: Mass spectrum of the neutral Higgs bosons, taking /ii = 200 GeV and A = 0.1, for: A = 200 GeV 
(left), A = 300 GeV (center), A ^ 0.1 GeV (right). 
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FIG. 11: Same as in Figure [TOl but taking A = 0.5. 



The piece of Lagrangian containing the fermion coupHngs of the charged Higgs bosons is given by: 



1 



V2VD 



C/3 



+ 



C/3 



+ 9i ]75 

S/3 



d + h.c, 



(18) 



where ((/>]^ )*, 02 are related to the physical charged Higgs boson states {H^ , H2 , H^) as follows: 



Hi 



(19) 



The UjkS denote the elements of the mixing-matrix that relates the physical charged Higgs 
bosons {H^ , H2 , H^) and the Goldstone boson G"*" (which gives mass to the W~^) with the fields 
(/)^, (l)^*, ^^'^ Ci* follows: 

f Un Ui2 Ui3 Uu \ / G+ \ 

U21 U22 U23 U2A 
U31 C/32 U34, 
\ Uii U42 f/43 U44 J 



it 

\ iV j 



(20) 



Then, the couplings udH^ , vilH^ are given by: 



9h+v,i 



i 



VD 



V2 



(^r+^r 75), qh, 

^'(1 + 75), 9h-u,i 



ud. 



ud 



VD 



V2 



VD 



V2 



^'(1-75) 



(21) 
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FIG. 12: Same as in Figure [TOl but taking A = 1. 
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FIG. 13: Mass spectrum of the charged Higgs bosons, taking A — GeV and A ~ 0.5, for: fii — 200 GeV 
(left), ^il = 400 GeV (center), /ii = 700 GeV (right). 
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FIG. 14: Mass spectrum of the neutral Higgs bosons, taking A — GeV and A = 0.5, for: = 200 GeV 
(left), /ii = 400 GeV (center), /ii = 700 GeV (right). In this scenario (B) only one massive pseudoscalar 
appears, as one can see from the limit A — (hence Bi — —A — 0) an extra Goldstone boson should appear, 
as discussed in (lOj . 



where and VJ" are defined as: 



a: 



ud 



mdtfs 



U2. 



i+l 



S/3 



niu cot/3 



ud 



rudtp h niu cot/3 



S(3 
^2,1+1 
S0 



Ul,i+1 
C/3 

C/3 



(22) 



One can see that the formulae in Eq. (10) become the couplings udHf, I'llHf of the MSSM when 
we replace U2,i+i — > sp and C/i,i+i — > — C/3[;2i]. The vertex udHf induces at tree-level the decay 
t — > b, which will be studied in the next section. 

B. The decay t H+ h 



In order to study this top quark BR we must consider both the decays t b for i = 1,2, 

because both modes could be kinematically allowed for several parameter configurations within 
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our model. The decay width of these modes takes the following form: 



^2 



,1 - + gh) 2 ^ — ~ ^^b- 



where A is the usual kinematic factor A(a, 6, c) = (a — 6 — c)^ — 46c and qi^ H+ = rn'^ f^+fm 



(23) 



Furthermore, we shall neglect the decay width for the light fermion generations. As we men- 
tioned before, if one replaces C/2,i+i — > and — > — C/3, the formulae of the decay width also 
reduce to the MSSM case: see, e.g., [28j]. In general the decay width for t — > b depends on 
the Superpotential parameters through the elements f^(i,2),j+i- When we choose an MSSM+ICHT 
scenario where spontaneous EWSB is dominated by the effects of the Higgs doublets, the decay 
width is practically similar to the one obtained in the MSSM. If we do not, results can be very 
different in the two models. We shall discuss this in the forthcoming numerical analysis. 



C. Numerical results for the decay t ^ H+ b in the MSSM+ICHT 

We explore several theoretically allowed regions of our MSSM+ICHT scenario and constrain 
these by using experimental bounds on the BR(t b). In the so-called "tauonic Higgs model" 



21I ]. the decay mode {H^ v^) dominates the charged Higgs boson decay width, and BR(t 

b) is constrained to be less than 0.4 at 95 % C.L. [21]. However, if no assumption is made on 
the charged Higgs boson decay, BR(t — > b) is constrained to be less than 0.91 at 95 % C.L. 



2l[. Conversely, the combined LEP data excluded a charged Higgs boson with mass less than 79.3 
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GeV at 95 % C.L., a limit valid for an arbitrary BR(i^"^ — > r"*" i^r) 

Thus, in order to conclude in this regard, we need to discuss all the charged Higgs boson 
decays following the steps of our previous paper [l(3]. In the present work, we shall evaluate all 
charged Higgs boson decays relevant masses below that of the top quark, thus including the modes 
T~^Vr,cs,cb,W~^Hf,W~^A^. In what follows, we want to find out whether a light charged Higgs 
boson (with mj{± < m^y±) is still allowed phenomenologically. As usual, we refer to our two 
benchmark scenarios. 

Scenario A. Remember that this scenario was defined by taking Bi = = 0, B2 = —A, and 
^2 = 100 GeV while for A we considered the values A = 0.1, 0.5, 1.0. In Figure [15] we present 
plots of: a) the BR(t bH^) vs. tan/?; b) the tan/? — plane; c) the BR(t — > bHi^) vs. 

tan/? and d) the tan/? - 771^^+ plane; for the case A = 0.1, taking A = 200, 300, 400 GeV. We 
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can observe that both modes t — > bH^ for i = 1,2 are kinematically ahowed. Also, we see that 
a charged Higgs boson with mass in the range 80 GeV < < 82 GeV and for 1 < tan /3 < 15 

satisfies the constraint BR(t bH^) < 0.4. Furthermore, from the plots of Figure [16] we can see 
that in this scenario the dominant decay mode is into t^i't, therefore they fall within the scope 
of the tauonic Higgs model, so that BR(t b) < 0.4 applies. This also happens for a heavier 

charged Higgs boson with mass in the range 125 GeV < Tnjj+ < 160 GeV and for 35 < tan (3 < 100. 
However, although this scenario is consistent with current experimental bounds from Tevatron on 
BR(t — > bH~^), it is excluded after one considers the LEP2 limits on the neutral Higgs sector via 



In Figure [T7] we present similar plots for the case A = 0.5, taking A = 200, 300, 400 GeV. We 
can observe that the mode t — > bH^ satisfies the constraint BR(t bH^) < 0.4 in the ranges 80 
GeV < 1^^+ < 115 GeV and 30 < tan/3 < 100. Then, from the plots of Figure [HI we see that 
in this range the dominant decay mode is into t'^i't, therefore they also fall within the realm of 
the tauonic Higgs model, so that BR(t — > H'^ b) < 0.4. Besides, this scenario is very interesting 
because we can have mjj± ~ myi/± and a light neutral Higgs boson, which could be consistent 
with the experimental LEP2 limits, as discussed in section II and shown in Table [H In Figure [19] 
we present the corresponding plots for the case A = 1.0, taking again A = 200, 300, 400 GeV. We 
can observe that in this case only the mode t bHf is kinematically allowed and satisfies the 
constraint BR(t bH+) < 0.4 in the ranges 140 GeV < < 160 GeV and 60 < tan/3 < 100. 

Scenario B. Recall that this was defined by taking B2 = ^2 = 0, Bi = —A, and /ii = 200 
GeV while for A we considered the values A = 0.1, 0.5. In Figure [20] we present plots of: a) the 
BR(t bH^) vs. tan/3; b) the tan/3 — plane; for the case A = 0.1, taking A = 200, 300, 0.1 
GeV. We can observe that the mode t bH^ is kinematically allowed for the case ^ = 0.1 GeV 
and, if we combine the results of Table IIIIl one can see that for charged Higgs boson masses in 
the range 84 GeV < fn^+ < 89 GeV and for 12 < tan /3 < 50 the model fulfills the constraint 
BR(t bH^) < 0.4. Again we have a charged Higgs mass ~ m\^±. Similarly, in Figure [2T] we 
present similar plots but for the case A = 0.5, taking A = 200, 300, 0.1 GeV. We can observe here 
that the mode t — > bHf' is kinematically permitted again for the case A = 0.1 and the model 
is fulfilling the constraint BR(t bH+) < 0.4 in the ranges 121 GeV < < 129 GeV and 

20 < tan/3 < 80. In Figure [22] we present only a plot of the tan(3 — m^+ plane, for the case A = 1.0, 
taking A = 200, 300, 0.1 GeV. We can deduce from here that the mode t — > bH^ is kinematically 
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FIG. 15: It is plotted: a) the BR(t bH^) vs. tan/3 (top-left), b) the tan/3 — plane (top-right), 

c) the BR(t bH2) vs. tan/3 (bottom- left), d) the tan/3 — "^^1^+ plane (bottom-right), in Scenario A by 
taking A = 0.1, for: A = 200 GeV (soHd), A = 300 GeV (dashes), A = 400 GeV (dots). 

forbidden, because mjj+ ~ 200 GeV in the range 1 < tan/3 < 100. Finally, in Figure [23] we present 
plots of: a) the BR(t bH^) vs. tan/3; b) the tanP — m^+ plane, for the case A = 0.5 and A = 0, 
taking fii = 200, 400, 700 GeV. We note here that the mode t — > bH^ is kinematically allowed 
and we see that a charged Higgs boson with mass in the range 121 GeV < Tn^+ < 131 GeV and 
for 1 < tan/3 < 70 satisfies the constraint BR(t bH^) < 0.4. Combining the results of Table 
IIV|, one can observe that only the small region 1 < tan /5 < 6 is allowed by the LEP collaborations' 
results. 

D. Decays of charged Higgs bosons in the MSSM+ICHT 

Let us now discuss the decay modes of the charged Higgs bosons within our model, which have 
an interest independently of whether these states are themselves produced in top decays. As usual, 
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FIG. 16: The figure shows the branching ratios of (top) and H2 (bottom) decaying into the principal 
modes in Scenario A, with A = 0.1 and /is = 100 GeV, for: A = 200 GeV (left), A = 300 GeV (center), A = 
400 GeV (right). The lines correspond to: (1) BR(ffi+ ^ r+Vr), (2) BR(iJi+ ^ cs), (3) BR(i7i+ -> cb), (6) 
BR{H+ T+Vr), (7) BR{H+ ^ cs), (8) BR{H+ ^ c6),(9) BR(iJ2+ W+H^), (10) BR(i72^ -> W+A\). 

we refer to our two customary benchmark scenarios. 

Scenario A. In Figure [TBI we present the BR's of the channels t^i^t, cs, cb, W~^Hi, W'^A'l 

for i = 1,2 as a function of tan/3 in the range 1 < tan/3 < 100 for the case A = 0.1, taking 
A = 200,300,400 GeV. When t bH+ is kinematically allowed for both i = 1,2, the dominant 
decay of the charged Higgs bosons is via the mode r+i/,-, with BR(//j^ t^Vt) ~ 1- We can 
observe that for H2 the decay mode M^+Aj* is dominant for tan {5 < 30, although the decay 
t — > is not kinematically allowed. In Figure [18] we present similar plots for , but now with 
A = 0.5 and again in this case the dominant decay mode is into r^Vr for the range tan (3 > 15. 
Then, from Figure [18] one gets that BR(ffj^ — > T~^Vr) ~ 1 when t bH^ is kinematically allowed. 
Similarly, in Figure[M]we present the corresponding plots for in the case A = 1.0. For A = 200 
GeV the dominant decay of the considered charged Higgs boson is the mode t'^i't, except in 
the range 15 < tan/3 < 35, where the decay channel W~^H^ is also relevant. For A = 300,400 
GeV the dominant decay of the charged Higgs state is via the mode t^Vt when tan/3 > 35, but 
for 12 < tan /? < 35 the decay channel W^H^ becomes the leading one whereas for the range 



25 




FIG. 17: It is plotted: the BR{t bH^) vs. tan/3 (left), the tan/? — plane (right), in Scenario A 

taking A = 0.5, for: A = 200 GeV (sohd), A = 300 GeV (dashes), A = 400 GeV (dots). 

2 < tan/3 < 9 the mode 1^"'"^]' is dominant. 

Now we discuss the decay modes of for the case A = 0.1. We can see in Figure [25] for A = 200 
GeV that the mode tb is dominant when tan /? < 2, but for 2 < tan /? the mode T^'*"^]' is the 
leading one. For the case A = 300 GeV there are three dominant decay modes: W'^H2 in the 
range tan /? < 2, W^Hf for 2 < tan /3 < 60 and tb when 60 < tan (3. For A = 400 GeV the relevant 
decay channels are: W^H2 in the range tan/3 < 2, W~^Hf for the two ranges 2 < tan/3 < 10 
and 20 < tan/3 < 40, W'^A^ for 10 < tan/3 < 20, tb when 40 < tan/3. In Figure [26] we present 
the corresponding plots for the BR's of the channels T^Ur,tb,W^Hj,W~^Al,W~^Z^ for 

z = 2, 3, J = 1, 2 as a function of tan/? in the range 1 < tan/3 < 100 for the case A = 0.5, taking 
A = 200, 300,400 GeV. When A = 200 GeV the dominant decay modes for H:^ are: W+A^ in the 
range tan/? < 25 and tb for 25 < tan/3. We obtain similar results for , except for tan/3 < 2, 
where the mode W^Hf is dominant. For the cases A = 300, 400 GeV the relevant decay modes 
of are: when tan/3 < 4, W^A\ in the range 4 < tan/3 < 25 and tb for 25 < tan/3. 

As for H2 , we obtain the same dominant modes as for in the range 4 < tan/3, because the 
decay channel W^H^ is relevant for tan/3 < 4. In Figure [271 we present plots for the BR's of the 
channels i?+ ^ r+i/^, tb, I^+if°, VF+A?, for i = 2, 3, j = 1, 2 as a function of tan (3 in the 
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FIG. 18: The figure shows the branching ratios of decaying into the principal modes in Scenario A, 
with A = 0.5 and /is = 100 GeV, for: A = 200 GeV (left), A = 300 GeV (center), A = 400 GeV (right). The 
lines correspond to: (1) BR(iJ+ r+Vr), (2) BR(i7+ ^ cs), (3) BR(iJ+ -> c5), (4) BR(i?+ W+Hl), 
(5) BR(i/+ ^ T4^+A?). 




FIG. 19: Same as Figure [T71 but taking A = 1. 
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FIG. 20: It is plotted: BR(i bH^) vs. tan/3 (left), the tan/? — plane (right), in Scenario B, taking 
A = 0.1 for: A = 200 GeV (solid), A = 300 GeV (dashes), A = 0.1 GeV (dots). 
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FIG. 21: Same as Figure but taking A = 0.5. 
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FIG. 22: The tan/3 - plane, in Scenario B, taking A = 1.0 foi-.A = 200 GeV (solid), A = 300 GeV 

(dashes), A = 0.1 GeV (dots). 

range 1 < tan/3 < 100 for the case A = 1, taking A = 200,300,400 GeV. As in the case A = 0.5, 
when A = 200 GeV the dominant decay modes for H2 are: in the range tan /? < 25, tb for 

25 < tan /?. We obtain similar results for , except for tan (3 <2 where the mode tb is dominant. 
For the cases A = 300, 400 GeV the relevant decay modes of H2 are: W~^Z^ when tan/3 < 4, 
W^jd]* in the range 4 < tan /3 < 25 and tb for 25 < tan /3. We get the following dominant modes 
for H^: W^H^ for 1.5 < tan/3 < 2, W^A^ when 2 < tan/3 < 25 while in two ranges tan/3 < 1.5 
and 25 < tan/3 the mode tb becomes the most important one. 

Scenario B. In Figure [251 we present the BK{H^ T'^Ur, cs,cb,W^ H^,W^ A^)^s as a function 
of tan/3 in the range 1 < tan/3 < 100 for the case A = 0.1 GeV, taking A = 0.1,0.5. When 
the mode t bHf is kinematically allowed, the dominant decay of the charged Higgs bosons 
is via the mode t^Vj- and it is obtained that BR(//j^ t^^t) ~ 1- We can observe that for 
A = 0.5 the mode PF^jdj* is dominant for 2 < tan /3 < 15. In Figure [29] we present BR(ifj^ — > 
r^i^T, cs, cb, W~^Hf, W^A^Ys as a function of tan /3 in the range 1 < tan /3 < 100 for the case ^ = 
GeV, A = 0.5, taking fii = 200, 400, 700 GeV. In these cases the decay t —>■ bH^ is kinematically 
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FIG. 23: It is plotted: BR(t bH^) vs. tan/3 (left), the tan/3 — plane (right), in Scenario B, taking 
A = 0.5 and A = for: fii = 200 GeV (solid), fii = 400 GeV (dashes), m = 700 (dots). 
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FIG. 24: Same as Figure [18] but taking A = 1. 
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FIG. 25: The figure shows the branching ratios of decaying into the principal modes in Scenario A, 
taking A = 0.1 and ^2 = 100 GeV for: A = 200 GeV (left), A = 300 GeV (center), A = 400 GeV (right). 
The lines correspond to: (11) BR{H+ th), (12) BR{H+ W+H^), (13) BR(iJ3+ ^ W+H^), (14) 
BR{H+ W+Hl), (15) BR{H+ W+A^l), (16) BR{H+ W+Z°). 

allowed and the dominant decay of the charged Higgs boson is via t^l^t for tan/3 > 20, since 
BR(//j^ t^Vt) ~ 1- In contrast, when tan/? < 20 the mode W^H^ is dominant. 

In Figure [30] we present the BR's of the relevant decay channels of the charged Higgs for 
A = 200, 300 GeV, taking A = 0.1. When tan/3 < 6 the dominant mode is tb. When 6 < tan/3 
though, the dominant mode becomes W^Z^. In Figure [3T] we observe that for the same previous 
values for A, but with A = 0.5, the dominant decays channels are: tb for tan /? < 40 and W^Z^ for 
tan (3 > 40. 

In Figure [32] we present plots for the BR's of the channels tb,W^ ,W~^ Al,W^ Z^ for 

i = 2, 3, J = 1, 2 as a function of tan/3 in the range 1 < tan/? < 100 for the case A = 0.1, 
taking A = 0.1,200,300 GeV. When A = 0.1 GeV the dominant decay modes for H2 are: tb for 
tan /? < 2 and W~^Z^ when 2 < tan /?. For the mode W^Z^ is the dominant one in the range 
1 < tan/3 < 100. For A = 200, 300 GeV the relevant decay modes of H2 are: tb when tan (3 <2 and 
W~^Z^ when 2 < tan/3. We obtain similar results for , but now the mode W~^Z^ is dominant 
when 6 < tan /3, and tb is the relevant mode for tan /3 < 6. In Figure [33] we present plots for the 
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FIG. 26: The figure shows the branching ratios of H2 (top) and (bottom) decaying into the principal 
modes in Scenario A, taking A = 0.5 and ^2 = 100 GcV for A = 200 GeV (left), A = 300 GeV (center), 
A = 400 GeV (right). The lines correspond to: (1) BR{H} r+Vr), (2) BR(iJ2+ ^ th), (3) BR(iJ2+ ^ 
W+H^), (4) BR(iJ2+ ^ W+A\), (5) BR{H+ W+Z"), (6) BR{H+ tb),(7) BR{H+ ^ W+H°), (8) 
BR{H+ W+H^), (9) BR(iJ+ W+Al), (10) BR(i73+ ^ W+Z"). 

BR's of channels tb, W+ ,W+ A'(,W+ Z° foi i = 2,3, j = 1,2 as a. function of tan p in the 

range 1 < tan/3 < 100, taking A = 0.5 for A = 0.1, 200, 300 GeV. When A = 0.1 GeV the dominant 
decay modes for H2 become: W^H^ in the range tan/3 < 1.5, tb for 1.5 < tan/3 and W~^Z^ when 
2 < tan/3. For the mode W+Z^ is dominant for ah tan/3. For A = 200, 300 GeV the relevant 
decay modes of H2 are: W^H2 when 2 < tan f3 and W^Z^ for 2 < tan /?. Similar dominant modes 
are obtained for , but now when tan/3 < 30 the mode W^Z^ becomes the principal one, and 
for 30 < tan/3 the mode tb is the relevant one. In Figure [Ml we present plots for the BR's of the 
channels T^Vr,tb, W^H^, W^A\, W~^Z° as a function of tan (3 in the range 1 < tan /3 < 100 

for A = 1, taking A = 0.1, 200, 300 GeV. For A = 0.1 GeV the dominant decay modes for are: 
W+Al in the range 2 < tan/3 < 40, tb for tan/3 < 2 and 40 < tan/?. For A = 200 (300) GeV 
the relevant decay modes of H2 are: tb when tan/? < 70(50) and VF"^yl5 when 70(50) < tan/3. In 
Figure [35] we present plots for the BR's of channels tb,W+ ,W+ Al,W+ Z^ for i = 2,3, 

J = 1, 2 as a function of tan /3 in the range 1 < tan /? < 100, taking A = 1 for ^4 = 0.1, 200, 300 GeV. 
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FIG. 28: The figure shows the branching ratios of decaying into the principal modes in Scenario B 
considering the cases: A = 0.1 (left), A = 0.5 (right), with A ^ 0.1 GeV and 200 GeV. The lines 

correspond to: (1) BR{H+ r+z/^), (2) BR{H+ cs), (3) BR{H+ cb), (4) BR{H+ W+H^), (5) 
BR{H+ ^14^+^0). 
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FIG. 29: The figure shows the branching ratios of decaying into the principal modes in Scenario B, 
with A = 0.5 and A = GeV for: ^ui = 200 GeV (left), /xi = 400 GeV (center), fj,i = 700 GeV (right). The 
lines correspond to: (1) BR{H^ r+Ur), (2) BR{H^ cs), (3) BR{H^ cb), (4) BR{Ht ^ W+H^). 




FIG. 30: The figure shows the branching ratios of decaying into the principal modes in Scenario B 
taking A = 0.1 GeV and /xi = 200 GeV, for: A = 200 GeV (left), A = 300 GeV (right). The lines correspond 

to: (1) BR{H+ T+Vr), (2) BR(ir+ ^ th), (3) BR(iJ+ ^ W+H^), (4) BR(F+ ^ W+Z°). 
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FIG. 31: Same as Figure [30l but taking A = 0.5. 

As in the case A = 0.5, when ^ = 0.1 GeV the dominant decay modes for H2 are: tb for tan (3 <2 
and W'^Z^ when 2 < tan/?. For the mode W'^Z^ is dominant for the entire tan/3 range. For 
the case A = 200 (300) GeV the relevant decay modes of are: W+H^ when 2 < tan /?, W+Z° 
for 2 < tan/? < 60(50), tb when 60(50) < tan/3 < 100(70). The dominant modes for the are: 
W~^Z^ in the range tan /3 < 60, and for 60 < tan (3 the mode tb is relevant. In Figure [551 we present 
the corresponding plots for the BR's of the channels — > tb, W~^Hj, W~^Ai, W'^Z^ for i = 2,3, 
J = 1, 2 as a function of tan /? in the range 1 < tan /? < 100, taking A = 0.5 for ^ui = 200, 400, 700 
GeV with A = GeV. The dominant decay modes for H2 are: W^H2 in the range tan/? < 1.5, 
tb for 1.5 < tan/3 < 2 and W^Z^ when 2 < tan/?. For the mode W^Z^ is dominant for all 
values of tan /3. For the case /Ui = 400, 700 GeV the relevant decay modes of H2 are: W~^Hi 
in the range tan/3 < 1.5, W+H^ when 1.5 < tan/3 < 2, for 2 < tan /3 < 60. Finally, the 

dominant mode for the is W^Z^. 

IV. DIRECT CHARGED HIGGS PRODUCTION AT THE LHC IN THE MSSM+ICHT 

We have found that, in some of the MSSM+ICHT scenarios envisaged here, light charged Higgs 
bosons could exist that have not been excluded by current experimental bounds, chiefly from LEP2 
and Tevatron. Their discovery potential should therefore be studied in view of the upcoming LHC 
and we shall then turn our attention now to presenting the corresponding hadro-production cross 
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A = 0.1GeV A = 200GeV A = 300 GeV 
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FIG. 32: The figure shows the branching ratios of ff^ (top) and (bottom) decaying into the principal 
modes in Scenario B, with A = 0.1 and = 200 GcV, for: A = 0.1 GeV (left), A = 200 GeV (center), 
A = 300 GeV (right). The lines correspond to: (1) BR{H+ tb), (2) BR{H+ W+H^), (3) BR{H+ ^ 
W+H^), (4) BR{H+ W+Al), (5) BR{H+ W+Z°), (6) BR{H+ tb),(7) BR{H+ W+H^), (8) 
BR{H+ W+H°), (9) BR{H+ W+A'^l), (10) BR{H+ W+Z^). 

sections via direct channels, i.e., other than as secondary products in (anti)top quark decays. 
As dealt with so far, if the charged Higgs boson mass rnr,± satisfies rnr,± < mt — mi,, where mt 

i i 

is the top quark mass and rrih the bottom quark mass, Hf^ particles could be produced in the decay 
of on-shell (i.e., Tt — > 0) top (anti-)quarks t bH^, and the charge conjugated (c.c.) process, the 
latter being in turn produced in pairs via qq annihilation and gg fusion. We denote such a 
production channel as qq, gg ^ tt ^ ^^^^ + c.c. (i.e., if due to (anti-)top decays) whilst we use 
the notation qq, gg — > tbH~ + c.c. to signify when further production diagrams are included^. In 
fact, owing to the large top decay width {Tt ~ 1.5 GeV) and due to the additional dia gram s which 



do not proceed via direct tt production but yield the same final state tbH- + c.c. 
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32 



331, 



^ Altogether, they represent the full gauge invariant set of Feynman graphs pertaining to the 2 — > 3 body process 
with a tbH~ + c.c. final state: two for the case of qq annihilation and eight for gluon-gluon fusion, see, e.g., 
Eq. (1.1) of 
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A = 0.1GeV A = 200 GeV A = 300 GeV 
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FIG. 33: Same as Figure [32l but taking A — 0.5. 



charged Higgs bosons could also be produced at and beyond the kinematic top decay threshold. 
The importance of these effects in the so-called 'threshold' or 'transition' regionf rn^i ~ rru 



was emphasized in various Les Houches proceedings 
so that the calculations of Refs 



are now implemented in HERWIG 



3S 
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3d, 



36 



371, 



35| as well as in Refs. 
32| fbased on the appropriate oo, gg thHf description) 
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421 and PYTHIA 



43 



44l |. A comparison between 
the two generators was carried out in Ref. [361. ^^"^ realistic simulation of production 



with m 



nit, as can well be the case here, the use of either of these two implementations is of 



paramount importance. 

Here, we use HERWIG version 6.510 in default configuration, by onsetting the subprocess 
IPROC = 3839, wherein we have overwritten the default MSSM/2HDM couplings and masses with 
those pertaining to the MSSM+ICHT: see Eqs. (j2ip -(|22 p . The production cross sections are 
found in Figures [37H39] and Figures H0Hi2] for our usual A (^i = 0) and B (^2 = 0) Scenarios, 
respectively, for various different choices of A and A. (See also Figure [43l illustrating the cross 
section dependence on tan /? for A = 0.5, ^ = and fii = 200, 400, 700 GeV.) 

The pattern of the cross sections reflects the usual dependence of the state of the 
MSSM/2HDM, ~ (mj cot/3^ + tan/3^)/(m^ +m^), induced by the Yukawa couplings inside the 
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A = 0.1GeV A = 200 GeV A = 300 GeV 
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FIG. 34: The figure shows the branching ratios of decaying into the principal modes in Scenario B, 
taking A = 1 and /ii = 200 GeV for: A = 0.1 GeV (left), A = 200 GeV (center), A = 300 GeV (right). 
The lines correspond to the modes: (1) BR(iJ+ r+Vr), (2) BR(i7+ tb), (3) BR(iJ+ W+Hf), (4) 
BR{H+ W+A"i), (5) BR{H+ ^ W+Z°). 

tbH~ vertex for the case of the , with a minimum at tan/3 ~ 6, as seen in Figures [37H391 our 
Scenario A. (In the last plot, one may appreciate also some peculiar interference effects between the 
'7^) 99 ^ tt ^ ^^^r + diagrams and the remaining ones, modulated by the top width.) The 
same can be said for the two heavier MSSM+ICHT states, and H^, with the minima shifted 
to lower tan/? values, the more so the heavier the Higgs boson, signalling that are the Yukawa 
couplings of these last two particles those differing most from the MSSM/2HDM limit. In the case 
of Scenario B, Figures H0Vl43^ the aforementioned coupling induced dependence upon tan/3 is only 
seen for the Hf^ state for very small ^'s (A = 0.1 and /xi = 200 GeV or even for A = GeV and 
/ii arbitrary). In all other parameter configurations, the trends differ dramatically from the 
usual (nearly parabolic) dependence typical of the MSSM/2HDM limit. The same can be said for 
all setups chosen in the case of the and H^. These peculiar patterns can be understood by 
using the analytic expressions given in Ref. [3] for the mass spectrum of the charged Higgs bosons 
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FIG. 35: Same as Figure [31 but taking A = 1. 
in the approximation R <^ 1. In this approximation, we can obtain in Scenario B that: 

m^± ~ 2niA/ sin 2/3, 



mrT± ~ 



Xv/{2V2R){Asm2(3 + 2ni) ±l/2{g'^ - A^)v^ 



(24) 



In the range 5 ^ tan/3 ^ 100, sin 2/3 decreases from ~ 0.1 to ~ 0.01. This explains the shapes of 
the curves showed in the Figures |40VI421 The parameter A just determines where the curves for 

m^± begin to fall. On the other hand, the pattern (m^ cot tan (m^ is suppressed 

^-2 



by a factor of 10 for rnjj± when A and tan (3 are large. 



Altoget 
MSSM in 



ler, by comparing th e qq, qq — > tbH^ + c.c. cross sections herein with, e.g., those of the 



451] or the 2HDM in 
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46( 1 ■ it is clear that the MSSM+ICHT rates can be very large 
and thus the discovery potential in ATLAS and CMS can be substantial, particularly for a very 
light , which may pertain to our MSSM+ICHT but not the MSSM or general 2IIDM. However, 
it is only by combining the production rates of this section with the decay ones of the previous 
ones that actual event numbers at the LHC can be predicted. 
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FIG. 36: The figure shows the branching ratios of H2 (top) and (bottom) decaying into the principal 
modes in Scenario B, taking A 0.1 and A ^ GeV for: ^1 = 200 GeV (left), fii 400 GeV (center), 
^ll = 700 GeV (right). The hnes correspond to: (1) BR(iJ2+ ^ tb), (2) BR(iJ2+ ^ W+H^), (3) BR(iJ2+ ^ 
W+H^), (4) BR(iJ2+ W+A°^), (5) BR{H^ W+Z°), (6) BR(i73+ ^ tb),{7) BR{H^ W+Hf), (8) 
BR(i?3+ ^ W+H°), (9) BR(iJ+ ^ W+A°), (10) BR(i73+ ^ W+Z"). 

V. CHARGED HIGGS BOSON EVENT RATES AT THE LHC IN THE MSSM+ICHT 

To illustrate the type of charged Higgs signatures that have the potential to be detectable at 
the LHC, we show in Tables IVl and IVll a summary of results for masses, LHC cross sections (u's), 
BR's and event rates. We focus on those cases where the charged Higgs boson mass is above the 
threshold for t — > H~^b. Thus, for Scenario A, all the entries for u's and BR's in Table IVl correspond 
to the second charged Higgs boson H^, while for Scenario B, Table IVTl shows the corresponding 
results for H^. All of these rates correspond to fi2 = 100 GeV for Scenario A and ^1 = 200 GeV 
for Scenario B. The cases where t — > H^b is allowed have been discussed previously in section HL 
We shall also assume an integrated luminosity of 10^^ pb^^. 

To illustrate these results, let us comment one case within each scenario. From Table [Vj we can 
see that for Scenario A, with A = 0.5, A = 200 GeV and tan /3 = 20, we have that H2 is heavier than 
mt — rrn,, with a mass rnjT+ = 390 GeV, this precluding top decay contributions, so that in this case 
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A = 200 GeV A = 300 GeV A = 400 GeV 
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FIG. 37: The figure shows the cross sections of H^2 3 the LHC through the channel qq, gg — > tbH + 
c.c. in Scenario A with A = 0.1 and for: A = 200,300,400 GeV, respectively. 

a{pp tbH+) ^ 1.2 X 10-2 pb, while the dominant decays are H+ tb, W+A^, W+Z'^, W+H^ 
which give a number of events of 180, 900, 12, 92, respectively. In this case the most promising 
signal is — W^A^. However, when A = 1.0 we have that all event rates get decreased, partially 
because the masses get enhanced. For instance mfr+ becaomes 545 GeV, however the number of 
events at the LHC for i?^ W^A^ can still be substantial, as it is about 349. 

Then, for Scenario B, we have that is already above the threshold for t — > H^b. So, for 
the declared values of the relevant parameters, we obtain that the lightest charged Higgs boson 
mass is Tnjj+ = 906 and 921 GeV for A = 0.5 and A = 1.0, respectively. In such a case only the 
decay —>■ tb can reach significant numbers for the LHC. We obtain a number of events of 912 
and 854, respectively. The other decay that has a large BR is — t'^Ut, but in all cases the 
number of events is at most of order 20-30, which seems quite difficult to be detectable at the 
CERN machine. 

Thus, we conclude that signatures in Scenario A are more diverse than for Scenario B. However, 

in order to reach quantitative conclusions we need to perform a study of signal versus background, 
but this seems more tractable if we concentrate in each mode individually, rather than while taking 
the general view that was attempted in this paper. 
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FIG. 38: The figure shows the cross sections of -f?i^2,3 LHC through the channel qq, gg — > tbH + 

c.c. in Scenario A with A = 0.5 and for: A = 200, 300, 400 GeV, respectively. 
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FIG. 39: The figure shows the cross sections of -f?i^2,3 LHC through the channel qq, gg — > tbH + 

c.c. in Scenario A with A = 1.0 and for: A = 200,300,400 GeV, respectively. 
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FIG. 41: The figure shows the cross sections of 3 the LHC through the channel qq,gg tbH~ + 
c.c. in Scenario B with A = 0.5 and for: A = 200,300,0.1 GeV, respectively. 
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A = 200 GeV A = 300 GeV A = 0.1 GeV 
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FIG. 42: The figure shows the cross sections of -f?i^2,3 the LHC through the channel qq, gg — > tbH~ + 
c.c. in Scenario B with A = 1.0 and for: A = 200,300,0.1 GeV, respectively. 




tanp tanp tanp 

FIG. 43: The figure shows the cross sections of -f?i^2,3 LHC through the channel qq, gg — > tbH~ + 

c.c. in Scenario B with A = GeV, A = 0.5 and for: /xi = 200, 400, 700 GeV, respectively. 
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VI. CONCLUSIONS 

After describing the entire Higgs mass and coupling spectrum, we have studied the fermion- 
charged Higgs bosons vertices within the context of an extension of the Minimal Supersymmetric 
Standard Model with an additional Complex Higgs Triplet (MSSM+ICHT). Then we have analyzed 
the decay t — bH^ (i = 1,2,3.) in the framework of this model. We have found that in some 
plausible scenarios of the MSSM+ICHT the experimental data allow regions in the plane M^± — 

i 

tan/3 for the lightest charged Higgs boson Hf' (even for H2) that are forbidden in the case of 
the MSSM. For instance, in this model it is possible to have a charged Higgs boson with mass 
mH± mw± which is not excluded yet by any of the current data. Furthermore, there are other 
regions of parameters where the top quark can decay into two charged Higgs states, so that plenty 
of these states could be produced at the LHC in (anti)top quark decays. Whenever the mass of the 
charged Higgs boson is larger than than mt, there is scope to exploit direct Hf" hadro-production 
through the associate mode qq, gg tbH~ + c.c. 

On the one hand, the detection at the LHC of charged Higgs bosons in the regions of parameter 
space accessible in the MSSM (or a 2HDM) would not contradict the MSSM+ICHT hypothesis. 
On the other hand, the observation at such a machine of several charged Higgs bosons would 
correspond to a model with a more elaborate Higgs sector than the MSSM (or a 2HDM), such as 
the MSSM+ICHT. 
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TABLE V: Summary of LHC event rates for Scenario A with p.^ = 100 GeV for an integrated luminosity of 
10^ pb~^, for several different signatures. 



0.5 



0.5 



0.5 



1.0 



1.0 



1.0 



(A,tan0) 



(200,5) 



(200,20) 



(200,50) 



(200,5) 



(200,20) 



(200,50) 



in GeV 



(118,740,790) 



(114,390,470) 



(98,290,370) 



(191,1047,1087) 



(185,545,610) 



(153,400,450) 



a(pp — » H2tb) in pb 



1.6 X 10-8 



1.2 X 10-2 



8.7 X 10-1 



3.4 X 10-® 



4.5 X 10" 



3.6 X 10" 



Relevant BR's 



BR 
BR 
BR 

BR 



BR 

BR 
BR 
BR 



BR 
BR 
BR 
BR 



BR 
BR 
BR 
BR 



BR 
BR 
BR 
BR 



BR 
BR 
BR 
BR 



tb) 



5.8 X 10" 



W+Af^ « 9.8 X 10-1 
W+Z°^ » 1.2 X 10-2 
W+Hf^ 4.2 X 10-3 



HI 



tb^ « 1.5 



X 10-1 
» 7.5 X 10-1 
Si 1.0 X 10-2 
^ 7.7 X 10-2 



tit 



tbj f5i 8.4 
■ W+A°^ 



8.2 X 10-2 
X 10-1 
» 2.4 X 10-2 
5i 4.8 X 10-2 



tb^ Ri 3.0 
W+A'l) : 



4.6 X 10-'^ 

X 10--* 
5i 9.8 X 10-1 
=s 1.3 X 10-2 



tb^ « 1.1 

■ W+Af^ : 

■ W+ZO^ : 



X 10-1 

« 7.7 X 10-1 

« 1.1 X 10-2 

K. 1.0 X 10-2 



■ tb^ Oi 8.4 



5.0 X 10-2 

X 10-1 
:s 2.6 X 10-2 
« 7.9 X 10-2 



Nr. Events 



180 
900 

12 
92 



7134 
73080 
2088 
4176 



50 
349 
5 
4 



1800 
30240 
936 
2844 
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TABLE VI: 

of 10^ pb-i, 



Summary of LHC event rates for Scenario B with /xi = 200 GeV for an integrated luminosity 
for several different signatures. 



0.5 



0.5 



0.5 



1.0 



1.0 



1.0 



(A, tan 0) 



(200,5) 



(200,20) 



(200,50) 



(200,5) 



(200,20) 



(200,50) 



in GeV 



(473,1304,1305) 



(906,1223,1225) 



(1206,1207,1424) 



(497,1838,1850) 



(921,1724,1738) 



(1434,1699,1713) 



a(pp — » H^tb) in pb 



1.8 X 10-2 



9.6 X 10-3 



2.0 X 10-« 



1.6 X 10-2 



8.9 X 10-3 



6.2 X 10-3 



Relevant BR's 



BR 
BR 
BR 
BR 



BR 

BR 
BR 
BR 



BR 
BR 
BR 
BR 



BR 
BR 
BR 
BR 



BR 
BR 
BR 
BR 



BR 
BR 
BR 
BR 



tbj « 9.J 
W+H°) 



1.7 X 10-2 

X 10-1 
« 4.2 X 10-^ 
« 5.4 X 10-" 



Ht 
Ht 



tb^ Ri 9.5 

W+ZO) 



3.9 X 10-2 
X 10-1 
» 1.4 X 10-5 
Si 2.1 X 10-'' 



Ht 
Ht 
Ht 
Ht 



tbj ^ 1.3 
W+H°) 



3.8 X 10-3 
X 10-1 
s 8.3 X 10-1 
« 3.6 X 10-2 



Ht 
Ht 
Ht 
Ht 



tbj Ri 9.7 
W+Z") ; 



1.6 X 10-2 
X 10-1 
a 4.0 X 10-^ 
5i 1.1 X 10-2 



Ht 
Ht 
Ht 
Ht 



T+vt) « 

tbj as 9.6 

11''+//-" 



3.9 X 10-2 
X 10-1 
« 5.8 X 10-6 

1.8 X 10-5 



Ht 
Ht 
Ht 
Ht 



3.8 X 10-2 
X 10-1 
5i 1.4 X 10-5 
W+H°^ 5.8 X IQ-^ 



tbj R! 9.6 



Nr. Events 



31 
1764 

1 



37 
912 






26 
1552 

18 



35 
854 





24 
595 
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